Upper Limits on the Extragalactic Background Light from the 
Very High Energy Gamma- Ray Spectra of Blazars 



M. Schroedter 1 

Whipple Observatory, Harvard-Smithsonian Center for Astrophysics, P. 0. Box 6369, 

Amado, AZ 85645-0097 

mschroed@cf a . harvard . edu 

present address: Dept. of Physics and Astronomy, Iowa State University, Ames, IA 

50011-3160 

ABSTRACT 

The direct measurement of the extragalactic background light (EBL) is diffi- 
cult at optical to infrared wavelengths because of the strong foreground radiation 
originating in the Solar System. Very high energy (VHE, E>100 GeV) gamma 
rays interact with EBL photons of these wavelengths through pair production. 
In this work, the available VHE spectra from six blazars are used to place upper 
limits on the EBL. These blazars have been detected over a range of redshifts 
and a steepening of the spectral index is observed with increasing source distance. 
This can be interpreted as absorption by the EBL. In general, knowledge of the 
intrinsic source spectrum is necessary to determine the density of the intervening 
EBL. Motivated by the observed spectral steepening with redshift, upper limits 
on the EBL are derived by assuming that the intrinsic spectra of the six blazars 
are oc E~ lA . Upper limits are then placed on the EBL flux at discrete energies 
without assuming a specific spectral shape for the EBL. This is an advantage 
over other methods since the EBL spectrum is uncertain. 

Subject headings: BL Lacertae objects: general — gamma rays: observations - 
diffuse radiation — infrared: general 



1. Introduction 

The measurement of the Extragalactic Background Light (EBL) is important for VHE 
gamma-ray astronomy as well as for astronomers modeling star formation and galaxy evolu- 
tion. Second only in intensity to the Cosmic Microwave Background, the optical and infrared 
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(IR) EBL contains the imprint of galaxy evolution since the Big Bang. This includes the 
light produced during formation and re-processing of stars. Current measurements of the 
EBL are summarized in Fig. 1. The optical to near-IR emission, peaking in the 1 /zm region 
and extending to 20 /zm is due to direct star light, while molecular clouds and dust reprocess 
the optical light and emit in the mid to far IR region producing the second peak, A ~ 20-300 
yum. Hauser and Dwek (2001) comprehensively reviewed measurements and implications of 
the cosmic infrared background. 

The optical to far-infrared EBL is difficult to measure because it is dwarfed by the 
much brighter foregrounds caused by night-sky glow, diffuse dust in the Galaxy, and the 
zodiacal light caused by interplanetary dust (Leinert et al. 1998; Hauser and Dwek 2001). For 
example, emission by the zodiacal dusk peaks in the 25 /xm region, orders of magnitude above 
the low EBL density in this waveband. In the case of ground- or rocket-based observations, 
instrumental emission also plays a significant role. This is complicated by the fact that 
the only characteristic upon which a detection of the EBL can be based is that it has 
to be distributed isotropically. These difficulties have precluded ground- and rocket-based 
measurements from detecting the EBL at all (Hauser and Dwek 2001). 

The measured flux of VHE gamma rays is attenuated by pair production with optical/IR 
photons. On extragalactic distances, most of this absorption occurs by the EBL (Nikishov 
1962); interactions with stellar and Galactic material and optical/IR photons are negligible 
(Dwek and Slavin 1994). Thus, if one somehow knows the initial gamma-ray flux, VHE 
astronomy is in a unique position to place limits on the density of the intervening optical/IR 
component of the EBL. 

The organization of this paper is as follows: existing measurements and constraints on 
the EBL are presented in Sect. 2. Then, the brightest flare spectra from each of the six 
blazars are presented in Sect. 3 together with the apparent spectral steepening with redshift. 
This is followed by a short review of pair-production absorption in Sect. 4 and is illustrated 
for the particular case of an assumed monoenergetic EBL in Sect. 4.1. Upper limits on the 
EBL density are derived in Sect. 5 and conclusions are given in Sect. 6. 

2. Measurements and Constraints on the EBL 

Direct measurements are possible in the two windows of least foreground around 1 
/im and >100 /iin (Hauser and Dwek 2001). Recently, the Cosmic Background Explorer 
(COBE) satellite with its two instruments, the Diffuse Infrared Background Experiment 
(DIRBE) and the Far Infrared Spectrometer (FIRAS), has detected the EBL at 140 fxm 
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and 240 jttm, see Fig. I. The possible detections at 60 /im and 100 /im (Finkbeiner et al. 
2000) are viewed as too high and are controversial, requiring revised galaxy evolution models 
with larger dust content (Blain and Phillips 2002). The FIRAS measurement (Fixsen et al. 
1998) shows that the IR EBL can be characterized between 125 /zm and 2000 mm by a 
modified blackbody spectrum. The isotropic optical and near-IR emission detected with 
the Japanese IRTS satellite is considerably higher than integrated light from galaxies and 
theoretical predictions (Matsumoto 2000). 

Lower limits on the EBL density are placed by adding the flux per unit area received 
from all galaxies down to a given flux limit. As galaxies are only one source contributing to 
the EBL, these galaxy counts represent a lower limit on the total EBL (Franceschini et al. 
1991; Armand et al. 1994; Pozzetti et al. 1998). In the mid-IR region, where the foreground 
is particularly bright, Elbaz et al. (2002) were able to place a lower limit on the 15 /im EBL 
density. 

Upper limits can be placed on the EBL from direct measurements with minimal back- 
ground subtraction (Hauser et al. 1998). Also, upper limits can be derived from fluctuations 
in the measured light distribution, see Kashlinsky et al. (1996); Kashlinsky and Odenwald 
(2000). 

2.1. Upper Limits from Observations of VHE Blazars 

BL Lacertae (BL Lac) objects are one type of active galactic nuclei that have a jet of high 
energy particles aligned with our line of sight. The nonthermal emission from the jet is highly 
variable and the mechanism for production of the VHE gamma-ray peak is under debate. 
In all models the intrinsic VHE spectrum of the source is smooth and concave downwards; 
no physical mechanism has been proposed that would produce an exponential rise with 
energy or emission-like line features. The VHE spectrum measured on Earth is modified by 
interactions of gamma rays with the EBL. For example, a 1 TeV gamma ray reaches the 
threshold for electron-positron pair production if it collides with a 1 eV optical photon. If a 
pair is created, this gamma ray is lost and causes an attenuation of the measured spectrum 
at that energy. Thus, the EBL density could be inferred from the measured spectrum if 
the intrinsic spectrum of the blazar were known. In addition, because of the coarse energy 
resolution of this method, a measurement of the EBL spectrum is only possible if the shape 
of the EBL spectrum is known in advance. 

Three methods have been used so far to derive EBL limits: 
1. Using multiwavelength observations, the intrinsic VHE blazar spectrum is predicted 
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through interpolation or fitting of an emission model. The EBL spectral shape is 
usually fixed in advance with only the absolute normalization of the EBL flux left as 
a free parameter (Stecker et al. 1992; Stecker and de Jager 1993; de Jager et al. 1994; 
Dwek and Slavin 1994; Biller et al. 1995; Guy et al. 2000). 

2. The deviation, or lack thereof, of the VHE spectrum from a predefined intrinsic spec- 
trum, e.g. a power-law, is used to place upper limits on the EBL density with a given 
spectral shape (Funk et al. 1998; Stanev and Franceschini 1998; Biller et al. 1998; 
Vassiliev 2000). 

3. For a given EBL scenario the intrinsic source spectrum is reconstructed. EBL models 
are ruled out, or limits are placed on the density, if they give rise to an unphysical 
source spectrum: such as that it must not rise exponentially and that it should be 
consistent with the X-ray synchrotron peak (Aharonian et al. 1999; Guy et al. 2000; 
Dwek 2001; Renault et al. 2001; Krennrich and Dwek 2003; Dwek and Krennrich 2005). 

The method described in this paper differs from the above in that: (1) The intrinsic 
source spectrum of all detected VHE blazars is assumed to be identical to each other. A con- 
servative upper limit to the intrinsic spectrum is interpolated from the observed correlation 
between the spectral index of the source during a flare state and its redshift. (2) Upper limits 
to the EBL intensity are calculated only at single energy points where the VHE spectrum 
has been measured, avoiding the need to specify the EBL spectral shape. 

3. Comparison of the Brightest Flare Spectra 

VHE emission has been detected and confirmed from six BL Lac objects: Mrk 421 (Punch 
et al. 1992; Petry et al. 1996), Mrk 501 (Quinn et al. 1996; Bradbury et al. 1997), 1ES 2344+514 
(Catanese et al. 1998; Tluczykont et al. 2003), 1ES 1959+650 (Nishiyama et al. 1999; Holder 
et al. 2003; Aharonian et al. 2003b), PKS 2155-304 (Chadwick et al. 1999; Aharonian et al. 
2004), and H 1426+428 (Horan et al. 2002; Aharonian et al. 2002a). 

For all six blazars, hysteresis and the rarity of simultaneous X-ray/ VHE spectra make 
it difficult to establish a temporal correlation between X-rays and VHE gamma rays. The 
sources are similar in that the X-ray synchrotron peak ranges from 10 17 to almost 10 20 Hz 
with peak luminosities between 10 43 and 10 44 erg s -1 sr -1 (Krawczynski et al. 2004). The 
black hole mass estimates are also relatively similar: 10 8 -10 9 M . It is possible that the 
reason for this similarity lies in the fact that only the brightest VHE blazars are seen by 
the imaging atmospheric Cherenkov telescopes at this time and that there is actually a 
continuous population of blazar spectra (Costamante et al. 2001). 
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The brightest flare spectra measured from each of the six blazars are shown in Fig. 2; 
they are increasingly steep with source distance. The spectra for Mrk 421 (Krennrich et al. 
2001) and Mrk 501 (Samuelson et al. 1998) are best described by a power law with exponential 
cut-off (Krennrich et al. 2001). The cut-off, at energy of about 4 TeV, can be attributed 
to gamma-ray absorption by the extragalactic medium. For 1ES 1959+650, the gamma-ray 
'orphan' flare is used (Schroedter 2004), a TeV flare during which the X-ray spectrum was 
not particularly bright (Krawczynski et al. 2004). For 1ES 2344+514, the brightest flare 
recorded thus far was recorded in 1995 by the Whipple 10 m telescope; the spectrum has 
been derived in Schroedter (2004); Schroedter et al. (2005). The spectrum of PKS 2155- 
304 has recently been measured by the HESS collaboration with high precision (Aharonian 
et al. 2004). Though the light curve shows variability, the spectrum can be well fitted by a 
single power law over the entire energy range while the power law index did not change with 
flux level. VHE emission from the furthest object, H 1426, is generally very low (Falcone 
2004), so that no flare spectra have been measured. However, as VHE emission during flaring 
episodes is much higher than during the quiescent state, most of the photons in the spectrum 
by (Petry et al. 2002; Aharonian et al. 2003a) may be associated with flares. 

The power-law spectral index fitted to each AGN flare spectrum is shown in Fig. 3. 
Because a significant cutoff is present in the spectra of Mrk 421 and Mrk 501, the spectral 
index is measured at two energies, 1 and 2 TeV, to estimate the systematic error arising 
from this choice. In both cases, the increasingly steep spectral index with redshift was fitted 
with a straight line, resulting in a spectral index at z = of -1.7+0.1 with x 2 /(d.o.f.) = 
22.3/4 when the spectral index is measured at 1 TeV and index of 2.0+0.1 (x 2 / (d.o.f.)=5.7/4) 
when the index is measured at 2 TeV. The flare spectra, after taking into account the redshift 
dependence, appear similar and hence we make the zeroth order assumption that the same 
process is responsible for the VHE emission in all these AGN. 

In the interpretation that the cut-off in the energy spectrum of Mrk 421 and Mrk 501 
is due to absorption by the EBL, note that for a power law fit with exponential cut-off, the 
power law index ranges from -1.9 to -2.1, with a cut-off energy around 4 TeV consistent for 
both blazars and independent of flux state (Krennrich et al. 2002; Aharonian et al. 2002b). 

3.1. Spectral Steepening Predicted by an EBL Model 

The literature contains various EBL models with various merits; a large fraction has 
been combined into different EBL scenarios by Dwek and Krennrich (2005). All of the 
scenarios that do not imply unreasonable intrinsic VHE spectra and most of the recent EBL 
models, produce similar gamma-ray opacities. Here, one particular EBL model is used as a 
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reference, put forward by Primack et al. (2001) and referred to as the Primack model in the 
following. This model is compared in Fig. 4 with EBL measurements and with the allowed 
EBL scenarios of (Dwek and Krennrich 2005). Although Primack et al. has proposed a 
number of models and the one chosen here is known to have slight problems above 60 jum, 
it will be used for illustrative purposes only, with the exception of giving an estimate of the 
optical depth already present for low energy gamma rays in Sect. 5. 

The attenuation of the intrinsic spectrum in the Primack EBL model is shown in Fig. 5 
at the distances of the six blazars. The figure shows the spectral slope, measured at 1 TeV, 
becoming increasingly steep with redshift. The steepening of blazar spectra with redshift 
was predicted by Stecker (1999). In the case of PKS 2155, the absorbed spectrum could be 
fitted with an E~ 3 spectrum over the energy region from 0.3 TeV to 3 TeV, assuming an 
initial E~ 2 spectrum. 

A linear fit to the spectral index steepened by the Primack model, measured at 1 
and 2 TeV, vs. redshift, is shown in Fig. 6. To determine the intercept of the linear fit, 
i.e. the intrinsic source spectral index, we assume that all blazars have the same intrinsic 
spectrum and adjust the level until the x 2 difference between measured and predicted indexes 
is minimized. The spectral index at zero redshift predicted in this way is -1.62 and -2.2, when 
the power law fits are performed at 1 TeV and 2 TeV, respectively. Considering the large 
degree of uncertainty in choosing an energy at which to measure the spectral index, Fig. 5 
implies that the intrinsic spectral index is on the order of -2 to -1.8. 

It should be noted that the detailed shape of the spectral steepening with increasing 
redshift depends strongly on the shape of the EBL. For example, the seven possible EBL 
scenarios given in Dwek and Krennrich (2005) allow a range of slopes for spectral steepening 
with redshift, with the Primack model roughly representing the average of the allowed EBL 
scenarios. 

Fig. 5 shows that in the Primack EBL model, there is a leveling-off at high energies, 
not a simple power law with exponential cut-off as supposed above. The cut-off energy in 
this case would correspond to the initial downturn at ~ 0.6 TeV. At higher energies, the 
flux is so far diminished that measurements with VHE telescopes do not yield a statistically 
significant result. Note that if the optical depth is independent of energy, as is the case 
for vF u (X) oc A -1 , no change occurs at all in the spectrum with increasing redshift. In the 
Primack model this is almost the case between 2 TeV and 5 TeV, corresponding to the 2 - 
5 /iin EBL. 
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4. Attenuation due to Pair Production 

Consider a single VHE gamma ray of energy, E, coming from a distant extragalactic 
source. If the gamma ray makes a collision with another particle, it will not propagate in 
a straight line and is lost to the observer who is several Mpc away. Pair production is the 
most likely type of inelastic collision in extragalactic space as the photon density, though 
varying with wavelength, is much higher than the matter density. For pair production to 
occur, the total energy available must be greater than 2m e . For VHE gamma rays with 
energy between 100 GeV and 20 TeV, the low energy photon must be in the range from 
10 eV to 0.05 eV. The absorption of VHE gamma rays is described by a very simple rate 
equation •jTeV + lir e + + e~. The number of VHE gamma rays, N, with energy E changes 
per unit time, dt, as 

= -N(E)\^ e+e -(E), (1) 

where A is the reaction rate defined below. Eq. 1 is solved by expressing dt in terms of the 
distance dl traveled by the photon: dl/dt = c. For small redshift, z << 1, the relation 
dl/dz = c/H holds so that the number of gamma rays at redshift z is given by 

N(z, E) = N (E) e~ H ° lz A(£) . (2) 

This defines the optical depth 

r(z,E) = H 1 z X(E), (3) 

with H = 71 ± 4 km s~ x Mpc -1 (Bennett et al. 2003). This formalism can been extended 
to cosmological distances (Salamon and Stecker 1998; Biller 1995). The optical depth of the 
furthest detected VHE blazar H 1426 at z — 0.129, calculated according to Eq. 3 and with 
the cosmologically correct version differs by less than 10%. As this difference is less than the 
systematic uncertainty of VHE blazar spectra, Eq. 3 will be used. The measured flux, F m , 
is then given in terms of the emitted flux, F e , by e~ T ( Em ^ F e (E rn ), where E m is the measured 
energy. This ignores the inverse-square law distance dependence; it is of no concern as it 
effects only the absolute flux level and does not have an energy dependence. 

The momentum distribution averaged pair production rate, A, in units of [T -1 ], is 

X(E) = (av rel ) = J d 3 p f(p) a(V^) v rel (4) 

= J dQ J p 2 dp f(p) <7(v/i) v rel (5) 

dx(l-x) p 2 dp f(p) a(y/I), (6) 

i J^?^ 



E(l-x) 
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where df2 = d(cos9)d(f), cosO = x is the angle between the incoming particles, is the 
azimuthal angle, and ^Js is the total energy. The relative speed between the interacting 
particles is v re i = c(l — x). The comoving photon momentum distribution of the EBL is 
given by f(p) [E^^L^^sr^ 1 ]. Eq. 6 can also be written in the more familiar way of particle 
density per energy per volume by substituting n(e) = 4irp 2 f(p) and with the replacement 
p — > e. The spin-averaged pair production cross section, a, is given by (Breit and Wheeler 
1934). From Eq. 6 it can be seen that the optical depth is independent of the gamma-ray 
energy if the EBL photon density, n(e), is independent of e: i.e. \(E) oc n only. This 
corresponds to an energy density dn/de oc e _1 . On a vF v plot, such an EBL spectrum would 
fall as A^ 1 and is quite possible in the optical / near-IR portion. 

4.1. Upper Limits from the Assumption of a Monoenergetic EBL 

The sensitivity of the VHE gamma-ray spectrum to absorption by the EBL is best 
illustrated by the toy model of a monoenergetic EBL, i.e. isotropically distributed photons 
with only a single energy. The absorption probability per unit length is given by 



where the constant n(e) is the co-moving photon density |X~ 3 ]. Note that the only depen- 
dence on the EBL energy e is through v e . The absorption probability per unit length is plot- 
ted in Fig. 7 for various EBL energies ranging from IR to UV: the density is taken uniformly 
as 1 cm~ 3 . Though gamma-rays are absorbed most efficiently when Ee = 0.93 x 10 12 eV 2 , 
there is significant absorption of gamma rays with energies from half to four times as much. 

In general, knowledge of the optical depth is not sufficient to unambiguously determine 
both the shape and magnitude of the EBL density. If the EBL flux is known over a finite 
wavelength region, an infinite number of shapes are possible because of the limited energy 
resolution of this method. The energy resolution is limited by the width of the pair pro- 
duction cross-section and the isotropic EBL photon distribution, see Eq. 6. If one assumes 
a shape for the EBL spectrum with only the overall normalization left as a free parameter, 
it is possible to determine the best fit EBL flux through a x 2 minimization between the 
measured and modeled optical depths. However, this still relies on knowing the intrinsic 
blazar spectrum, which is not easily accomplished. This method has been widely used to 
derive upper limits for power-law EBL spectra. 

To avoid specifying the EBL spectral shape altogether, upper limits on the true EBL 
density are derived here by assuming the EBL to be monochromatic. Physically, we know 




(7) 



- 9- 



zH^ 1 



that the EBL spectrum is extended. If we make that assumption that the EBL is monochro- 
matic, the contribution of other wavelengths to the absorption of gamma rays is ignored. 
Hence, the EBL density at that single wavelength will have to be larger to reproduce the 
optical depth, than if other wavelengths were allowed to contribute to the absorption. 

Suppose we start with a known measurement of t(E) from a source at z — 0.031, 
illustrated in Fig. 8 (left) for three arbitrarily chosen values of t(E). Then, we want to 
know: What EBL density does it take to reproduce this amount of absorption? The answer 
to this question depends on what photon wavelength is performing the absorption. The 
least EBL density is required if it absorbs the gamma-ray most efficiently, that is, when 
Ee = 0.93 x 10 12 eV 2 . The monochromatic EBL density corresponding to this case can be 
expressed in terms of the measured optical depth, r m , using Eq. 3 

c f +1 

n(e)- J ^ dx(l-x) a(v e ) (8) 

TrnHp 

->n(e) = Tl z - , (9) 

f Ll dx (! - x ) a ( V e) 

so that the upper limit (UL) on the true density is given by 

.O.gSeVTeV, imiEUtu 
n UL ( ) = 14x ' 1Q % 9m2 . (10) 

Fig. 8 (middle) shows the EBL flux corresponding to this case and Fig. 8 (right) shows the 
corresponding optical depth produced by this single monochromatic EBL line. 

Consider then the case of many adjacent optical depth measurements, shown in Fig. 8 
(left) at only three energies. As before, the monochromatic EBL flux is calculated in turn 
from each optical depth (middle). The total absorption, the sum of the three optical depths, 
is larger than was measured at the three energies (right). Hence, the monochromatic EBL 
flux is an overestimate of the true value required to achieve the measured attenuation. 



5. EBL Upper Limits from Blazar Spectra 

Determination of the amount of absorption present in one measured spectrum requires a 
priori knowledge of the intrinsic source spectrum. Though the mechanism for VHE gamma- 
ray production has been modeled, which model to choose is still unknown. In addition, other 
important details are left open as well: What is the electron or proton spectrum? What is 
the strength of the magnetic field? What is the opacity in the vicinity of the source for 
gamma-rays to escape? 
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To place upper limits on the EBL density, the following simple assumption will be made: 
The intrinsic spectrum of all blazars are identical to a power law with index -1.8; this has 
been motivated in Sect. 3. If the VHE spectrum has not been measured down to 0.2 TeV, 
the flux at that energy is predicted from a power law fit extrapolated down from higher 
energies. The small amount of EBL absorption present already at 0.2 TeV is calculated from 
the Primack model. 

The assumed spectral index of -1.8 is very hard and cannot extend indefinitely. By 
using a power law fit from higher energies, where the measured spectrum is steeper due to 
absorption, to extrapolate down in energy to 0.2 TeV will likely overestimate the expected 
flux from the source. The assumed intrinsic spectral index is approximately implied by the 
physics of TeV blazars. The intrinsic spectra of TeV blazars are expected to peak in the 
multi-TeV region (Stecker et al. 1996); hence the intrinsic spectral index should be somewhat 
harder than -2. Detailed modeling show this to be the case for Mrk 421 and Mrk 501 
(Krawczynski et al. 2002; Konopelko et al. 2003). Dwek and Krennrich (2005) show that 
after accounting for absorption by the EBL, the intrinsic spectra, in a vF v representation, 
peak between 0.5-1.2 TeV for Mrk 421, between 0.8-2.5 TeV for Mrk 501, and between 1-5 
TeV for H 1426. 

The expected flux in the absence of EBL absorption from each blazar is shown in Fig. 9 
by a dashed line. The measured flux levels have been adjusted to the distance of Mrk 421 
using an inverse square law to show that the power output of the sources could be similar 
if the sources have identical geometries. The optical depth is the ratio of expected flux to 
measured flux with the addition of the optical depth already present at 0.2 TeV 

t = In (F e /F m ) + t(z, E = 0.2 TeV) (11) 

and this does not dependent on the overall normalization of the flux level. The optical depth 
calculated in this way for the six blazars is shown in Fig. 10. Also shown in the figure is 
the optical depth calculated in the Primack model for comparison with the measured values. 
A comparison of the measured optical depth for Mrk 421 with the allowed EBL scenarios 
of (Dwek and Krennrich 2005) is given in Fig. 11. Here, the Primack model lies roughly 
in the middle of the allowed EBL scenarios. The measured optical depths of Mrk 421 and 
Mrk 501 are steeper in the 1-5 TeV region than some of the EBL scenario and would favor 
a larger EBL density in the mid-IR region, as in the IR spectrum derived by Malkan and 
Stecker (2001). The measured optical depths of the other sources are flatter, showing that 
this sample of six blazars is not completely homogeneous. 

The monochromatic EBL flux derived from each blazar spectral energy point under the 
above assumptions is shown in Fig. 12. For comparison, the value of the monochromatic flux 
as calculated from the Primack model is shown in the figure as well. For the most part, the 
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monochromatic flux values derived here agrees with the Primack model; this is also true for 
all of the other reasonable EBL scenarios. 

Differences in the derived monochromatic flux are likely caused by the intrinsic flare 
spectra of the sources being different from each other and not being pure power laws. Because 
of the large distances to the blazars, it is highly unlikely that non-uniformities in the EBL 
give rise to the observed differences unless an IR galaxy is in the line of sight. The limits 
derived from Mrk 421 and Mrk 501 are the lowest. To bring them higher and in line with 
the other sources, their intrinsic spectra would need to be slightly harder than -1.8 in the 
low energy region. If the spectra were harder then curvature with a peak in the 1-2 TeV 
range would make the flux limits more consistent with those from the other sources. The 
systematic uncertainty in absolute energy calibration is 10%, directly corresponding to a 
10% uncertainty in the EBL wavelength in Fig. 12. 

Fig. 13 shows the derived 98% confidence level upper limits together with optical EBL 
measurements and limits shown previously in Fig. 1. The range of upper limits from the six 
sources is indicative of the systematic error of this method. The systematic error stemming 
from the optical depth at E — 0.2 TeV is not shown; it is at the level of 5% at 1 TeV and 
decreases with wavelength. The mid-IR limits on the EBL are monotonously decreasing with 
wavelength and above the limits inferred from fluctuation-analysis of the EBL by Hauser and 
Dwek (2001). The upper limits in the optical to near-IR are in conflict with the detections 
claimed by Matsumoto (2000) and also with the HHH scenario by Dwek and Krennrich 
(2005). An increase in the upper limits in this region would be achieved if the intrinsic source 
spectrum peaked in the 0.7-2 TeV region or if the optical depth at 0.2 //m is substantially 
higher. 

Matsumoto (2000) argues that his measurements are inconsistent with galaxy evolution 
models and much higher than what can be accounted for by the observation of galaxy 
populations. Already at 2.2 fim the flux is higher than that claimed by (Wright 2001) 
from COBE and 2MASS data. They found that zodiacal emission modeling as their main 
uncertainty and responsible for the large fluxes. 

The limits on the EBL derived here are compared with others, derived also with VHE 
gamma-ray observations, in Fig. 14. Other upper limits are generally lower; this is due to the 
assumption of a continuous EBL spectrum. In some cases, the upper limits are in conflict 
with measurements using optical methods. 
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6. Conclusion 

The EBL is difficult to measure directly because of bright foreground radiation originat- 
ing within the Solar System. VHE gamma-rays are very penetrating and shine through most 
of the dust in the Galaxy and the Solar system. On cosmological distances, pair production 
with optical and infrared radiation attenuates the gamma-ray flux. Therefore, if the VHE 
source spectrum is known, the EBL density can be inferred from the measured attenuation in 
the spectra. Unfortunately, the known extragalactic sources of VHE gamma rays are highly 
variable and modeling of the source spectrum is difficult. To estimate the intrinsic source 
spectrum, the strongest flare spectra of the six sources were compared. These were found 
to be similar once absorption by the EBL was accounted for. Therefore, the zeroth order 
assumption was made that the VHE flare spectra of BL Lacs are actually the same; with 
differences in attenuation arising solely due to the different distances to the objects. This 
simple, and almost experimental, way to infer the intrinsic spectrum at the chosen energy 
of 0.2 TeV, results in a spectral index expected from the physics of TeV blazars, see for 
example (Stecker et al. 1996; Krawczynski et al. 2002; Konopelko et al. 2003). 

The method developed here to derive upper limits on the EBL does not hinge on an 
assumed EBL spectral shape. In detail, it consists of the following three steps: (1) The optical 
depth for gamma-rays was derived at discrete spectral energy points from the difference 
between the measured gamma-ray spectrum and the assumed intrinsic spectrum. With 
better understanding of the intrinsic VHE blazar spectra, the optical depth, and hence the 
resulting upper limits, can be determined more accurately. (2) The density of EBL photon 
with a single wavelength was found that reproduces the measured optical depth at a given 
gamma-ray energy. The EBL wavelength was fixed so that the pair production is most 
efficient, requiring the smallest density of EBL photons. (3) The derived monochromatic 
EBL density represents an upper limit to the true EBL density. The true EBL spectrum 
is extended and absorption of a gamma ray is produced by EBL photons over a range of 
wavelengths, thereby requiring a lower density than a monoenergetic EBL to produce the 
same optical depth. 

The optical depths measured from the six blazars fall in line with predictions from a 
range of EBL models. In particular, for Mrk 421 the measured optical depths represent 
about the average of the allowed EBL scenarios given by Dwek and Krennrich (2005). The 
upper limits derived here are substantially above other limits that were derived assuming 
a continuous (or stepwise continuous) EBL, see Fig. 14. However, this does not mean that 
the limits derived here are weak; to properly compare them with EBL models, upper limits 
must be calculated from the EBL models first. This was done for one EBL model, shown in 
Fig. 4 to represent an average of the allowed EBL scenarios, and the limits derived from the 
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model and the measurements are in agreement. 

The constraining power of the upper limits depends on the EBL spectrum: if the spec- 
trum is steeply falling with wavelength, the upper limits are more conservative at shorter 
wavelength than if the EBL spectrum is flat. The unambiguous and statistically correct de- 
termination of 98% CI upper limits from combining the results of the six blazars is difficult 
because an average intrinsic blazar spectrum was used, derived from the observed spectral 
steepening with redshift. The upper limits derived here should be considered as a whole and 
not from one blazar alone. 

The derived upper limits question one particular set of EBL measurements (Matsumoto 
2000) that are very high in the near infrared waveband and which galaxy formation models 
are typically not able to reproduce either. 

In principle, one could extend the method described here to derive more constraining 
upper limits on the EBL density by reiterating a process of calculating upper limits and then 
using the upper-limit spectrum as an input to a second iteration of calculating upper limits. 
To do better and measure the EBL density, requires accurate knowledge of the intrinsic and 
the absorbed source spectra over at least two decades of energy. Reduced error bars on the 
measured spectrum will be achieved with the new generation of telescopes coming online 
now. However, the mechanism for production of VHE gamma rays, and hence the intrinsic 
source spectrum, is still under considerable debate. Unless a theoretical model is accepted or 
a "standard candle" at cosmological distances is identified with a well understood gamma- 
ray spectrum, the spectra from all these sources will only provide upper limits on the EBL 
flux. Only through further observation at all wavelengths will the gamma-ray production 
mechanism be understood. The GLAST satellite, to be launched in 2006, will cover the 
energy range from 0.1-100 GeV. Together with the low threshold of about 100 GeV for the 
new arrays of imaging Cherenkov telescope such as HESS, VERITAS, and CANGAROO, 
simultaneous flux measurements can be performed over almost six decades in energy. 

The author acknowledges the award of a Smithsonian Predoctoral Fellowship and the 
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Fig. 1. — Spectral energy distribution of the EBL. The top axis indicates the most likely 
energy of the partner photon to participate in pair-production with the EBL. The flux 
measurements and limits are taken from (Hauser et al. 1998) (open circles), (Finkbeiner et al. 

2000) (open squares). Upper limits from fluctuation analysis by Hauser and Dwek (2001) 
using the data from Kashlinsky et al. (1996); Kashlinsky and Odenwald (2000) (stars). 
Tentative detection by (Miville-Deschenes et al. 2002) (open diamonds). (Lagache et al. 
1999) (filled diamonds), (Dwek et al. 1998) (open triangle pointing down), (Fixsen et al. 
1998) (shaded region) near the CMB. Galaxy counts corresponding to lower limits are shown 
with open triangles pointing up and are at 0.2/im from (Armand et al. 1994), from 0.36 /zm 
to 2.2 fim from (Pozzetti et al. 1998), and at 15 fim from (Elbaz et al. 2002). In the infrared, 
the filled circles are from (Wright and Reese 2000), filled triangle pointing up from (Wright 

2001) , and filled triangles pointing down from (Cambresy et al. 2001). In the optical, the 
filled squares are from (Bernstein et al. 2002). The grey line in the optical/infrared is from 
(Matsumoto 2000). The cosmic microwave background (CMB) is shown by a dashed line. 
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Fig. 2. — The very high energy flare spectra of six AGN. The shaded region shows the fit 
and 68% CI for 2 independent parameters. Mrk 421 from Krennrich et al. (2002) (open 
diamonds), Mrk 501 from Samuelson et al. (1998) (filled diamonds), 1ES 2344+514 (filled 
circles), 1ES 1959+650 (open circles), PKS 2155-413 in 2003 (left triangles) (Aharonian 
et al. 2004), H 1426+428 from Petry et al. (2002) (filled squares) and Aharonian et al. 
(2003a) (open squares). For clarity, the flux level of 1ES 2344, 1ES 1959, and H 1426 has 
been reduced by a factor of 2, 10, and 10, respectively. 
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Fig. 3. — Spectral index of power law fits to the hardest measured flare spectra versus 
distance of the VHE blazars. The index was measured at 1 TeV (left) and at 2 TeV (right) 
for the curved spectra of Mrk 421 and Mrk 501. A linear fit to the data is shown in each 
case. 
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Fig. 4. — The Primack model with CMB flux added (dashed line). Left: with measurements 
and upper limits as shown in Fig.l (left graph). Right: with the allowed EBL scenarios given 
in (Dwek and Krennrich 2005) (thin lines). 
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Fig. 5. — Illustration of power law spectra measured at 1 TeV from the attenuation predicted 
by the Primack model. The attenuation, e~ T (solid lines), is shown for the blazars (top to 
bottom, solid lines) Mrk 421, Mrk 501, 1ES 2344, 1ES 1959, PKS 2155, and H 1426. The 
dashed lines show the power-law spectra at 1 TeV. These become steeper with increasing 
source distance. No attenuation occurs at z — 0, the horizontal line at the top of the figure, 
e~ T = 1, corresponds to the intrinsic source spectrum and could have any reasonable shape. 
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Fig. 6. — Spectral index of power law fits to the hardest measured flare spectra versus 
distance of the VHE blazars. The index was measured at 1 TeV (left) and at 2 TeV (right) 
for the curved spectra of Mrk 421 and Mrk 501. The prediction of the Primack model for the 
steepening of the spectra is shown by line with points. The predicted index was evaluated 
at each source distance and produces an almost linear relationship. For the Primack model, 
the spectral index at zero redshift is a free parameter and is derived from a best-fit to the 
data. The reduced \ 2 f° r the fit in the left figure is 5.7 and 4.3 for the right figure. 
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Fig. 7. — Absorption probability of a VHE gamma ray by a monoenergetic EBL of energy e 
= 0.01 eV (123 pim, solid), 0.1 eV (12.3 /im, dashed), 1.0 eV (1.23 /im, dash-dot-dot), and 
10 eV (123 nm, dotted). The EBL photon density is 1/cm 3 . 
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Fig. 8. — Illustration of the method to derive an upper limit on the EBL density from a mea- 
surement of the optical depth. From the measured optical depth (left), the monochromatic 
EBL intensity is calculated (middle). The (right) graph shows the optical depth produced 
by the monochromatic EBL line. Note that the optical depth derived from each monochro- 
matic line overestimates the measured optical depth everywhere except at the measured 
gamma-ray energy. The total absorption is the sum of the three lines. 
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Fig. 9. — The very high energy flare spectra of six AGN. The flux of each source is adjusted 
with an inverse square law to the same distance as Mrk 421. Shown by the dashed line in 
each case is the assumed source spectrum with differential index -1.8; they are normalized 
to the measured, or predicted, flux at 0.2 TeV. The top left graph indicates with an arrow 
the effect of the optical depth in attenuating the flux at two energies. 
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Fig. 10. — The optical depth (data points with error bars) derived according to Eq. 11 from 
the differences in Fig. 9 between measured and assumed flux. For comparison, the optical 
depth derived from the Primack model is shown by a solid line. 
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Fig. 11. — The optical depth (data points with error bars) derived according to Eq. 11 for 
Mrk 421. Superimposed are the optical depth derived from the Primack model (thick line) 
and optical depths derived from the allowed EBL scenarios given in Dwek and Krennrich 
(2005) thin lines. 
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Fig. 12. — The monoenergetic EBL flux derived from the spectra of AGN (points with error 
bars). Also shown for comparison is the monochromatic flux (thick line) derived from the 
Primack EBL model (thin dashed line). 
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Fig. 13. — 98% confidence level upper limits of the EBL flux (points with upper limit error). 
Measurements and limits from other publications are shown in gray, see also Fig. 1. The top 
axis indicates the most likely energy of the partner photon to participate in pair-production 
with the EBL. 
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Fig. 14. — The gray region shows the 98% confidence level upper limits of the EBL flux 
derived in this work. Upper limits from other publications are as follows: Stecker and de 
Jager (1993) (dash), Dwek and Slavin (1994) (line with stars), Biller et al. (1995) (thick line: 
dot- dot- dash), Funk et al. (1998) (thin line: dot- dot- dash), Stanev and Franceschini (1998) 
(thin line), Biller et al. (1998) (thick line), Guy et al. (2000) (dot-dash), Vassiliev (2000) 
(thick dotted line), Dwek (2001) (line and single point with squares, Renault et al. (2001) 
(line with diamonds) . The (dotted parallelogram) shows the EBL flux estimate of de Jager 
et al. (1994). All flux values have been rescaled to H = 71 km/s/Mpc. 



